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H Abstract. In this article, we will get the estimate of the expected distribution of critical 

C ~ 1 values of Gaussian SU(2) random polynomials as the degree large enough. The result is a 

direct application of the Kac-Rice formula. 

o 

Q 1. Introduction 

Random polynomials and random holomorphic functions are studied as ways to gain insight 
for problems arising in string theory and analytic number theory [TJ |H1 [TT] . In [T2] , Kac studied 

i and determined a formula for the expected distribution of zeros of some real Gaussian random 

polynomials. His work was generalized to complex random polynomials and random analytic 

p I functions throughout the years, we refer to [2j [3l [5], [Til HI] f° r more backgrounds and results. 

1.1. SU(2) polynomials. When the random polynomial is defined invariant with respect to 

S some group action, the problem can turn out to be particularly interesting, we refer §2.3 in 
[TT] for examples. In this article, we will study a special family: the Gaussian SU{2) random 
polynomial. This is of particular interest in the physics literature as the zeros describe a random 
p<| spin state for the Majorana representation (modulo phase) on the unit sphere [§]. 

J> Given a probability space and {a^JLi a collection of i.i.d complex Gaussian random 

variables with mean and variance 1 on it, the family of 577(2) random polynomials is defined 

oo as ' 

O (1) Pn(z) -- 

(N 

— Although this polynomial is defined on C, we may also view it as an analytic function on 

CP 1 =CUoo with a pole at oo. 

Various properties of the zeros of random SU (2) polynomials have been studied such as 
the distribution of zeros and the two points correlation function [§]. First, zeros of this 
polynomials are uniformly distributed on S 2 = CP 1 with respect to the Fubini-Study metric, 
i.e., the average distribution of zeros is invariant under the SU(2) action on CP 1 [TT]. To be 
more precise, let's denote 

zGCP 1 : p„(z)=0 

as the empirical measure of zeros of Gaussian SU(2) random polynomials and define the pairing, 

(Z Pn ,4>)= ^( z )> where 0e C*°° (CP 1 ) 

zeCP 1 : p„(z)=0 
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We define the expectation, 

(EZ Pn ,<t>) := E(Z Pn ,<l>) = jf^ \ I e ~^de ao ■ ■ ■ dl a „ 

where dl aj = ^da,j A ddj is the Legesgue measure on C. 
Then the expected density of zeros is [2], 

in the sense that, 

E(Z Pn ,4>) = n I (jxjJFS-, where 4> G C° 
Jcv 1 

where uips is the Fubini-Study metric on CP 1 [5]. 

We can also study the two points correlation function of zeros of SU(2) polynomials and its 
scaling property. We define the two points correlation function as [2], 

K n (z,w) :=E{Z Pn {z)®Z Pn {w)) 

such that for any smooth test function <j>i(z) (g> <p2(w), we have the pairing, 

(K n (z, w),Mz) ® foM) = E {(Z Pn , 0i)) {{Z Pn , (j) 2 )) 

If we scale the two points correlation function by a factor then we have, 

z to (sinht 2 + i 2 ) coshi — 2t sinhi ^ ^ 1 
\/n ^/n sinhi 3 

where £ = L— ^1 and |z — w\ is the geodesic distance of z and to on CP 1 . It's easy to see, 



2 

z w . 2 



* n (^,_) = t --*- + O (O-st->0 

which implies zeros repel each other. We refer to [2j |9] for more details. 

1.2. Main results. In this article, we will study the expected distribution of nonvanishing 
critical values of \p n \ as n large enough. 

Note that the modulus \p n \ is a subharmonic function, thus there is no local maxima; local 
minima are all zeros and thus nonvanishing critical values are obtained only at saddle points 
[B] . Hence the expected density of nonvanishing critical values of \p n | we study in this article is 
in fact the expected density of values of saddle points of \p n | . 

The nonvanishing critical values of \p n \ are obtained at points, 

(2) {z&C: p' n =0 and p n ^ 0}. 

For a random polynomial p n , it has no repeated zeros almost surely, which implies that the set 
([2| is almost surely equivalent to, 

(3) {zeC: 74 = 0}. 

i.e., (nonvanishing) \p n \ and p n have the same critical points almost surely. 

Hence, we will first get the expected density of critical values of p n in Theorem [T] as a direct 
consequence, we can apply the polar coordinate to get the expected density of nonvanishing 
critical values of \p n | in Theorem [5] 

We denote the empirical measure of critical values of p n as, 

( 4 ) c p„ = 

z: p' n (z)=0 
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We now define the pairing, 

(5) {C Pn ,4>)= E V0(x) G C c °°(K 2 ) 

z: p^(z)=0 

where C^°(R 2 ) is the space of smooth functions on IR 2 with compact support. 

We denote D Pre (x) as the expected density of critical values of p n in the sense that, 



(6) E(C Pn ,<f>) = / </)(x)I]> Pn (x)d£ x , V4>(x) € C™(R 2 ) 

Jc 

whereas dl x is the Lebesgue measure of C. 

Those definitions also apply to the empirical measure of the nonvanishing critical values of 

\Pn\, 

( 7 ) C \p,i\ = <w 

z: p'„=0 

which is a measure defined on the nonnegative real line M. + . 
We define its expectation as, 



(8) (EC| P „|,0) :=E(C K |,0) = / 4>(xp\ Pn \dx, V <f>(x) € C C °°(M+) 

Jo 

where dx is the Lebesgue measure on M. 

In this article, we will first get the exact formula for the expected density D Pri in the Propo- 
sition [I] by the Kac-Rice formula (see section then we study the asymptotic behavior of D Pn 
as n — > oo. Our main results are, 

Theorem 1. The expected density U> Pn of the empirical measure C Pn of the critical values of p n 
satisfies the estimate, 

l _ —\ x \ 2 i r 1 

(9) D Pri = + -/ e- {s - slo&s ^ 2 ds + o{l) as n^oo 



7T _ 

for x G C. 

It seems that the growth D Pri is as |x| is away from as n — > oo, below is the computer 
graphic of this function with the variable t := \x\, 



As proved in Proposition [TJ the density B> Pn d£ x only depends on i.e., the modulus of |f> n |, 
thus we can rewrite it as D Pii (|x|)|a;|(i|a;|d0 under the polar coordinate. If we integrate on 9 
variable, then 

Jo 

will be the density of critical values of \p n \. Thus as a direct consequence of Theorem [11 we get, 
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Theorem 2. The expected density Oi Pn i of the empirical measure Cu i 0/ £/ie nonvanishing 
critical values of \p n \ satisfies the estimate, 

(10) %„|(a!) = 2< ^— e —l +2x f e-( s - slogs)x2 ds + o(l) as n -> 00 

/or x € M+ . 

The growth of D| p „| is also of order 1/x as x is large, but there is a peak when x is small, the 
following is the computer graphic, 




10 15 20 



1.3. Further remarks. First note that our setting is different from the one in [4j. For example, 
in [4], critical points of SU(2) polynomials are defined to be the points 

{z £ CP 1 : V Pn = 0} 

where V = J| — jz^z is the smooth Chern connection on the line bundle 0(n) — > CP 1 with 
respect to the Fubini -Study metric and p n is a global holomorphic section of the line bundle 
0{n) — > CP 1 [5]. By choosing such smooth Chern connection, the expected distribution of 
critical points are also invariant under the SU(2) action [4.. But in this article, the critical 
points are defined by the usual derivative 

dp n 



{zt 



dz 



= 0}. 



In fact, the derivative ^ is a meromorphic flat Chern connection on 0(n) CP 1 with a pole 
at 00. Under this setting, the expected density of critical points is not SU(2) invariant, we refer 
to 10J for more details. 

Our second remark is as following. In [7], the authors studied the expected density of non- 
vanishing critical values of the pointwise norm of Gaussian random holomorphic sections of the 
positive holomorphic line bundle over compact Kahler manifolds. Now let's briefly explain the 
main result in [7] and compare it with Theorem [2] Take Gaussian ST/ (2) random polynomials 
(sections) p n for example. We equip the line bundle 0(n) —> CP 1 with a Hermitian metric 
where <fi — log(l + |z| 2 ) is the Kahler potential of Fubini-Study metric. Then the 



h n = 



-ntf> 



pointwise /i-norm of the holomorphic section \p n | ^ 
hence the critical points of |p n |/i» is defined as, 

E„ ={zeQ 



\Pr, 



2 is global defined on 



and 



d\p 



'n\h n 



= 0} 



We define the (normalized) empirical measure of critical value of |jj n |/i n as , 



-|p„| h « 



"Ip»I*» 

\zgS„ 
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which is also a measure defined on R + . 

Then the expectation of C| p „k„ satisfies the estimate, 



(11) W.C Mhn =x [2x^-4 + 8e--)e- x +0(-), x£ 



n 



as n large enough. In fact, this estimate is universal: it holds on any Riemannian surfaces [7J. 

Thus the (normalized) density EC| Pn | h „ is exponent decaying as x large enough which is quite 
different from the behavior of (non- normalized) density EC| Pll | in Theorem [2j This is mainly 
because of the connection we choose: the usual derivative 4- in this article is a meromorphic flat 



connection on CP 1 with a pole at oo while in [7J, the proof of (11 ) relies on a choice of smooth 



Chern connection V' = 



d nzdz 

i+M 2 



Acknowledgements: The first named author would like to thank S. Zelditch for suggesting 
this problem and many helpful discussions. 



2. Kac-Rice formula 

In this section, we first review the Kac-Rice formula for a stochastic process, referring to 
[U H21 H3] for more details. Then we generalize the formula to the expected distribution of 
critical values of p n . 

The Kac-Rice formula is as follows: let f(z) be a real valued stochastic process indexed by a 
compact interval I cK. Then the Kac-Rice formula for the expected number of zeros is, 

E#{ze/: f(z) = 0} = ( ( \y\p z (0,y)dydz 



where p z (0, y) is the joint density p z (x, y) of (/, /') evaluated at (0, y). If / is a Gaussian process, 
then the joint density p z (x, y) is determined by the covariance matrix of (/, /') [T]. 

The proof of this formula is explained in more details in £Q . The idea of the proof is based 
on the following observation, 



\dz 



#{zel: f(z) = 0} = js (f{z))\f(z)\c 
We take expectation on both sides to get, 

E#{z E I : f{z) = 0} = / / S (x)p z (x,y)\y\dxdydz 



\y\Pz(0,y)dydz 
Thus the expected density of zeros of / is given by, 



(12) E 



I E M = ( f \y\Pz(o,y)dy) dz 

W:/(*)=0 J J 



If f(z) is a complex stochastic process indexed by a compact complex domain, then the above 
formula reads, 



(13) El ]T 6 A = (J c \y\ 2 p z (0,y)d£ v ^jd£ z 



.ZEI: f(z)=0 
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where d£ y and di z are Lebesgue measures on C. Compared with ( 12 1, we get \y\ 2 since a 
1-dimensional complex random process is a 2-dimensional real random process. In fact, this 
formula is based on the definition of the delta function, 

#{z € I : f{z) = 0} = J 5 (f(z))~df A df = J S (f(z))\f\ 2 dl z 

The formula arises when we take expectation on both sides. 

2.1. Kac-Rice formula: Revisit. In this subsection, let's get the formula for the expected 
density of critical values of a (real or complex) stochastic process / by the method of Kac-Rice. 

For simplicity, let's first consider a smooth real stochastic process / £ C°°{I) where I is a 
compact subset in R. 

Let 8 c R be a compact subset. Let's denote the set of critical values in Q as, 

C e = {zel: /Wee, /'(*)= 0}. 
Let's denote the measure [i(x)dx on as, 



VzGCe / 



in the sense that, 

(j}/j.(x)dx 



where <f> is any smooth test function defined on O. 
Then we have the following lemma, 

Lemma 1. Let's denote p z (x, y, £) as the joint probability of (/, /', /"). The 

fi(x)dx = I / / 0, £)dl;dz I dx 



where dx, g?£ and dz are Lebesgue measures on 
Proof. First note that, 



/ee, /'=o /ee, /'=o Jl 



By taking expectation on both sides, 



\fe&, /'=o / 



Pz(x,y,i)x{xeQ}4>{x)&{y)\S,\dydidzdx 

Pz (x,0,0\^\dCdz\ <j){x)dx 
4>(x)fi(x)dx 

□ 
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In the proof of Lemma [T] we have assumed I and being compact subsets in HL But the 
proof of Lemma[I]can be generalized to the SU(2) random polynomials p n which are a collection 
of complex Gaussian stochastic processes indexed by C. 

The generalization of to be C only requires picking up a sequence of discs centered at the 
origin with radius m € {1, 2, . . .} and taking limit in weak sense. And the generalization from 
I to C is the same. 

However, we do need to modify the pairing by choosing the test functions <f>{z) in the smooth 
compact supported space (E 2 ) in order to change the order of the integration on C. Following 
the proof of Lemma [T] we have, 

Lemma 2. The expected density of critical values of p n is, 

(14) ® Pn di x =([ [ \i\ 2 Pz {x, 0,0d£ c de z )d£ x 




where d£ x ,d£^ and dl z are Lebesgue measures on C and 

(15) ^ X ^ = ^^ z ^[-((°\^ 

is the joint density of {PmP'mPn) where A z is the covariance matrix of (p n ,p' n ,p'^). 

The proof of this formula is the same as the one in Lemma [l] Note that we get |£| 2 in 
the formula since p n is a 1-dimensional complex Gaussian process which is a 2-dimensional 
real process. Moreover, p„ is a Gaussian process, hence the joint density p z (x,y,£) is uniquely 
determined by the covariance matrix of {p n ,p' n ,p'n), we refer to [TKHIl] for more details. 

3. Proof of main Theorems 

3.1. The Density U> Pn . In this subsection, we will derive the exact formula for D Pri based on 
Lemma [2] We prove, 



(16) « r , = ^l ^^e-^^r 



Proposition 1. The expected density of the empirical measure of C Pn is given by the formula, 

n - 1 
7T Ji r" 
Thus U> Pn is a function only depending on \x\. 

Proof. By Lemma [2j in order to compute the expected density of critical values of p n , we first 
need to compute the covariance matrix of (p n , p' n , p'JJ . 

By definition, the covariance matrix of the Gaussian process (PmP'mPn) is given by [US], 

(Efa&d WnP^) E(pW 

A = E(p4) E(pX) E(pX) 
\E(p„K) Eip'rji") E(M 

The covariance kernel for the Gaussian process p n is, 

E(p n {z)p^{w)) := n„(z, to) = (1 + zw) n 

Then we can express each entry in the covariance matrix as, 

dU n (z,w) 



E{p n p n ) = n„(z, z), E(p' n p n ) = 



d 



z 



E , „ x _ d 2 U n (z lW ) _ d 2 U n (z,w) 
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WnPn) 



Some straight forward computations show that the covariance matrix is, 



/ 1 



A z = (i + |z| 2 r 

Hence, 
(17) 

Now we denote, 



nz 

i+M 2 



n(n — l)z 2 



i+M 2 (i+M 2 ) 2 (i+kl 2 ) 3 

.t(w-l)z 2 2n(ra-l)z+(n-l)n 2 z|z| 2 2n(ra-l)+4ra(n-l) 2 \z\ 2 +n 2 (rt-1) 2 |z| 4 
V(l+|z| 2 ) 2 (l+|z| 2 ) 3 (1 + M 2 ) 4 



det A z = (1 + \z 



2\3n 



2n 3 - 2n 2 



Then by direct computations, we rewrite, 
(l + |z| 2 ) 2 " //V 



det A. 



(1 + |2| 2 ) 6 

(n 3 -n 2 )z 2 
(1 + M 2 ) 4 



(1 + N 2 ) 4 
re 

(i+M 2 ) 2 



We expand this expression and further rewrite Q z (x, £) as, 



(18) 



1 



2(l + |z| 2 )™ 



Vn 2 - nz 2 x + . I g(l + |z| 2 ) 2 



— n 



+ 2(n|z| 2 + l)|a:| 5 



By Lemma [2j the expected density of critical values of p n is given by the formula, 



(19) 



oJ x ) 



1 



;-Q»(*,f) 

det A 2 



|£| 2 <^4- 



(20) 



Let's integrate £ variable first. Plug (|18|) into (|19|), we can rewrite (19 1 as, 



1 



7T° 



e (i+m 2 > 



^N 2 



det A, 



where A z is the following integral in £ variable, 

1 



K z = exp 



A', 



2(l + |z| 2 )« 
—I 

—n 

(n 2 — n) 2 



n 2 — nz x + 



1 



l£l 2 <% 



Change variables £ -> ^J,^ ! + M 2 ) 2 to get, 



exp 



2(1 + \z\ 2 ) n 



(i + kl 2 ) 8 7c 

Further change variable, £ — > \Jn 2 — nxz 2 + £ to get, 
(n 2 -n) 2 /" f |£| 2 



\J n 2 — nxz 2 + £ 



A z = 



(i + ki 2 ) s y c exp i 2(1+1^) 



(if, 



which equals to 



AT, = 



(n 2 — n) 2 



exp 



ICI 2 
2(1 + kl 2 )' 



|> £ — \/r? — nx 
y\e + (n 2 -n)\xz 2 \ 2 )d£ 
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Integrate it out, we have, 



K z = [2(1 + M 2 rV - n)M a |z|* + 4(1 + \ z \r n ] 



Now we change variable r = 1 + \z\ , rewrite 



and 



Kz = ^ n ) 2 [2r n (n 2 - n)\x\ 2 (r - l) 2 + 4r 2 



det A z = r 3n - 6 (2n 3 ~ 2n 2 ), e " C?H^ |x|2 = g-^^^M 2 



Plug these two lines back into the formula of (20) and use the polar coordinate di z — \drd0, 
integrate on 9 variable, we can rewrite D Pii as, 

(21) D P = — f ^-^"(^l) 2 N 2 + 2r 2 " e _^ N2dr 

y ' tt Ji r 3n + 2 

There are two parts in the numerator, we use integration by part to simplify the first part. Note 
that 

de -^±iM 2 = e-^^^M 2 ^-"- V " n)(r - l)\x\ 2 ]dr, 
then the first part is equal to, 

f0 ° (n 2 ~ n)r n (r ~ l) 2 \x\ 2 e _ ^-i) + i lxl 2 ^ 

j.3n+2 



n 




1 




7T 




n 




1 




7T 




n 




1 


7T 



r 

oo 



n+1 



Tl n+1, rc(r-l) + l i ,2 



Hence the density (21) is further simplified to be, 



n — 1 

which completes the proof. 

3.2. Proof of Theorem [l] Now we turn to the proof of our main Theorem [l] 
We denote t= \ and 

Vn{t) = n(r ~ r l ) + 1 = nt-- 1 - (n - 

then we have £ € [0, 1] and y n (t) & [0, 1] with y n (0) — and y n (l) = 1. 
Substitute n< " r ~^ +1 by t/„(t), we rewrite D Pn in Proposition (HI) as, 

^ n-1 /-"^(t) 



7T 

n — 1 

7T 



L e -»-(t)i*r dr 



2/ n (t)e-^WI a l 2 df 



where in the last step, we change variable t -. 

3 calculate 

9n{\x\ 2 ):= /V'-WN'di, 



The trick to estimate D Pii is to calculate 



□ 
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By integration by part, we have, 

g n {\x\ 2 ) = /Ve-^WM'dt 
Jo 

=e-^ 2 + f 1 ty' n (t)\x\ 2 e-y^W 2 dt 
Jo 

ri 

=e-^ 2 +n(n~ l)\x\ 2 / (t™" 1 - t n )e~ v " ^ dt 
Jo 

=e-W 2 +n\A 2 (nt"- 1 - (n - l)t n )e- y ^ x \' dt - n\x\ 2 / /" '< " l " ' "«// 
=e-l^ 2 + n|x| 2 / y^e""^!^* - |x| 2 /i„(|x| 2 ) 



N 3 + ^^ Dpn _ N X(N 2 ), 



=e 

where we denote 

(22) MM 2 ) := n / r^-w-WH di. 



Thus, 



(23) D p „ = tfnU U|a + ( * 2 

717T \ |X| Z 

We claim 

lim 5„(|a;| 2 ) = 1. 

n—¥oo 

This is quite straight forward, as Ve <E (0, 1), we rewrite, 



g n (\ x \ 2 ) = / e~y^ 2 dt 



Since y n (t) — > uniformly on [0, 1 — e] as n — > oo, thus 



lim / e -*»WM*di= / lim e-^WI^tft = 1 - e 



For the second integration, since y n (£) > on [0, 1], we have f, e ''"•WW 2 dt < e Hence we 

l-e< lim / e - y '^ x \ 2 dt< lim / e^C'I'I'dt^l 

As e is chosen arbitrarily, letting e — > + yields the claim. 
Now we estimate ( 23 ) to be, 

n - 1 ( 1 - e-l^l 2 X 

= 

(24) 



p„ I | | 2 + MM ) + 0(1) 

mr \ xr 



l-e-l^l 1 
= — ^^ + "MM 2 ) + o(l). 

as n — > oo. 

We now turn to estimate h n (\x\ 2 ). Change variable s = t n , h n will be rewritten as 
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where 

It's easy to check that 



z n (s) = —ns n + (n — l)s. 
z n (s) < Z n +\{s) 



for any fixed s G [0, 1]. 

Thus we have z n (s) monotone increasing to — (s — slogs) as n — > oo, hence, /i„(|x| 2 ) will 
satisfy 



lim hj\x\ 2 ) 



-(s-s log s)|a| ^ s 



This will give us the estimate 



h n (\x\ 2 ) = / e -( s - sloss ^ 2 ds + o(l) 
Jo 



as n — > oo. 



Hence we further estimate (24 1 to be, 



1 - e-W 1 



tt\x\ 2 tt J 
which completes the proof of Theorem [T] 



+ - / e -(>->l>>*')\*\ ds + o(l) asii4 



[l 
[3: 

[4: 
[5: 

[7. 

[»: 

[9 
[10 
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